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Abstract 

A new method based on the Monte-Carlo calculation on the lattice is proposed to 
study the Casimir effect in the compact lattice U(l) theory with Wilson action. We 
have studied the standard Casimir problem with two parallel plane surfaces (mirrors) 
Q\ ', and oblique boundary conditions on those as a test of our method. Physically, this 

boundary conditions may appear in the problem of modelling of the thin material 
films interaction and are generated by additional Chern-Simons boundary term. 
This approach for the boundary condition generation is very suitable for the lattice 
formulation of the Casimir problem due to gauge invariance. This method can be 
simply generalized on the case of more complicated geometries and other gauge 
groups. 
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1 Introduction and general motivation 

o : 

\Q , During the last few years Casimir effect has attracted much attention due to the great exper- 
imental and theoretical progress in studying of this phenomenon. This macroscopic quantum 
effect plays crucial role in nanophysics, micro-mechanics, quantum optics, condensed matter 
physics, material science, and also it is very important for different models of boundary states 
. in hadron physics, heavy-ion collisions, and cosmology. 

Nowadays, many theoretical methods for calculation of Casimir effect are proposed. Un- 
fortunately, most of these methods are based on the fixed boundary conditions or external po- 
tentials. Such a simplification can lead to problems in gauge invariance, renormalizability and 
locality. Moreover, there are some problems with different thermal corrections to Casimir force: 
different methods of calculations of the Casimir effect predict different corrections. Another 
part of the problem is calculation of Casimir force for complicated forms of boundary surfaces. 
Typically, various approximate methods (like the proximity force approximation method pQ, 
[2]) are used in the case of curved surfaces, but it is still unclear now whether those approx- 
imations are correct and for what tasks they can be applied. And finally, analytic methods 
for Casimir effect calculation are very complicated and strongly dependent on the shape of 
surfaces. Practically, Casimir effect was studied analytically only for cases of plane, spherical 
and cylindrical surface forms. More complicated tasks have not been studied very well by this 
moment but such cases typically appear in experiments. 

Another very interesting and important problem is the calculation of Casimir effect for other 
gauge groups (for example, SU(2) and SU(3)) and for fermion fields; both are very important 
for particle physics and for cosmology. 

Based on arguments discussed above, it seems a very important task now to create a general 
method for calculation of Casimir effects which would work well for different shapes of boundary 
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surfaces, for different fields, at non-zero temperature and density and under other external 
factors. It means that such a method should be formulated very generally for working in 
different coupling regimes and different external conditions. And we believe that only direct 
lattice calculations in quantum field theory can meet all these requirements. 

In our paper we consider the simplest Casimir problem with two parallel planes as a test of 
our method. The existence of the analytical answer for this problem is the additional motivation 
for such a choice of bound surfaces. This analytical answer assists us in fitting procedure for 
our numerical results. 

A crucial obstacle on the way of the realization of any Casimir problem on the lattice is the 
following. If the Casimir energy is a reaction of the vacuum on the presence of the boundary, 
what is "the boundary" in terms of lattice formalism? In other words, what is an observable 
quantity corresponding to such boundary? In fact, the answer is non-trivial and we devote the 
first part of the paper to this question. In the second part, we discuss the lattice algorithms 
and numerical results. And the last part is a conclusion. 

2 Chern-Simons boundary conditions and Casimir effect 

Casimir effect is a reaction of the vacuum on boundary condition. The spectrum of vacuum 
fluctuations depends on the boundary conditions. Changing of the boundary conditions leads 
to changing of the spectrum of vacuum fluctuations and so to generating of the corresponding 
Casimir force on the boundary. In the standard quantum field theory formalism, such changing 
of spectrum of vacuum fluctuations can be described, for example, by means of Green function 
method pQ. This approach is a very powerful tool for studying many essential Casimir tasks 
Pp. Unfortunately, the application of this analytical method to the case of more complicated 
shape of the boundary surfaces is not so easy due to calculation difficulties. Our aim is the 
creation of the numerical method for the Casimir effect calculation directly from the quantum 
field theory action. The lattice formalism looks very attractive for this role but manifestly we 
can not base in our approach on the separation of vacuum modes corresponding with boundary 
from the full spectrum of vacuum fluctuation. In lattice formalism we work in Euclidian space 
and deal with full spectrum of vacuum fluctuations and can not easily snatch out vacuum fluc- 
tuations corresponding to some boundary conditions. We need some very delicate approach for 
separation of vacuum fluctuation modes that preserve gauge invariance of our lattice formalism. 
Fortunately, such an approach to Casimir problem was proposed recently 

This approach is based on a very elegant idea coming from some unique properties of the 
Chern-Simons action in three dimensions [5j[6]. Let us consider electro-magnetic fields in 3+1 
dimentions with the Maxwell action and additional Chern-Simons action given on 3-dimentional 
integral on the boundary surface S: 

S = ~\J d ' x V" ~ ^<fd 3 se^ p n a A,(x)F up (x), (1) 

where e a ^ v P is the Levi-Civita tensor and n a is the normal vector to the boundary surface S, A 
is a real parameter. 

Let us consider now the simplest form of boundary surface S, namely two parallel infi- 
nite planes placed at the distance R from each other. The Chern-Simons formulation of this 
canonical Casimir problem was studied analytically in series of works [TOl ITT] . We will use this 
analytical answer for the fitting of our numerical data. 

In the case of plane form of the boundary surface S the Chern-Simons action in ([TJ has the 
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Figure 1: /(A). 



following form: 

S cs = \j (6(x 3 ) - 5(x 3 - R))e 3 ^A^x)F up (x)d 4 x. (2) 

where, in our formulation of this Casimir problem, normal vectors to the planes are turned 
in the opposite directions. This choice of the normal vector orientation corresponds to our 
renormalization procedure based on the connection between open and closed Casimir problems. 

If the parameter A is small, electro-magnetic fields obviously don't feel any boundary and 
are free. What happens if the parameter A becomes large and tends to infinity and fields 
dynamics on the boundary surface S is determined by Chern-Simons action? Let us consider 
the equation of motion obtained from the action ([T]): 

DA" + X(5(x 3 ) - 5(x 3 - R))e^ up A a d v A p = 0. (3) 

At A — > oo, it is easy to obtain from ([3]) a corresponding boundary conditions on the surface S: 

E n \ s = 0, #|||s = 0, (4) 

where E n and H\\ are normal and longitudinal components of electric and magnetic fields 
correspondingly. These conditions mean the nulling of the energy flux of the electromagnetic 
field through the surface. 

This phenomenon corresponds to the well-known property of the Chern-Simons theory in 3- 
dimension, namely to topological vortex (strings) generation [7JE]. As it was shown by Witten 
[9], Chern-Simons theory, which is a quantum field theory in three dimensions, is exactly 
solvable by using nonperturbative methods and the topological vortices play a leading role in 
this process. It was shown that the partition function of this theory depends on the topology 
of the surface S and gauge group only, and tends to the zero if A tends to infinity. There are no 
propagation modes in this case and if the boundary surface S is closed, the gauge fields inside 
and outside this surface become separated from each other by the topological vortex fields on 
the surface. In the case of the finite value of A, we have a non-trivial interaction between 
topological vortices and electro-magnetic waves [8]. For our practical purposes it is enough to 
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know that the parameter A plays the role of the effective regulator of the penetrability of our 
boundary surface for electro-magnetic wave. 

The analytical answer of Casimir energy per unit area for two planes is given by Green 
function method is the following [Till [TT] : 



7T 2 



E Cas = "72oF/(A)> (5) 



where function /(A) : lim A ^ 00 /(A) = 1 (Fig. 1) can be written as: 



7T 4 \X 2 + 1 

here the polylogarithm function Li±(x) is defined as 

00 „& 



Li 4 (x) = V" ^— = — / k 2 ln(l - xe~ k )dk. 
k=i k 2J o 



In our paper we will study the Casimir energy per unit area behavior at small A, which is 
the following: 

fi*. = - S ^ + 0(V) (6) 

3 Wilson "bag" and numerical lattice simulation of Maxwell- 
Chern-Simons theory 

In our paper we use the four- dimensional hyper-cubical lattice and the simplest form of the 
action for the U(l) gauge fields (so-called "Wilson action"): 

X flU 

where the link and plaquette variables are defined as: 

TT _ igA^a 

Up^ x ,fiu — Ui yXtfl Ui )X +ji tl/ u^ x+ p^ u\ x v . 

Here a is the step of the lattice and the lattice parameter (3 = 1/g 2 . Physical quantities are 
calculated in the lattice formalism by means of field configuration averaging, where the field 
configurations (the set of all link variables) are generated with the statistical weight e~ Sw . 

We have clarified in previous sections that additional the Chern-Simons action describes 
Casimir effect. In order to find a lattice description of Casimir interaction between boundary 
surfaces let us consider Wilson loop, which describes the interaction of charged particles. Wilson 
loop can be written in QED as 

ig § A^dx^ . . 4 

W c = e c = e i!J»A»dx\ ( 7 ) 

The exponent in ([7j) is the additional term to the action. This term describes the interaction 

of the field with the current J^(x) = g § S(x — £)<i£ M of charged particle. Configuration 

c 
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Figure 2: Wilson bag for two plane surfaces. 



averaging of Wilson loop (W(R,T)) (where R and T are dimensions of the loop) converges in 
the limit T — > oo to: 

(W(R,T))^Ce- v ^ T , 

where V(R) is the energy of interaction between charged particles. The same method can be 
used for calculation of Casimir energy by means of Chern-Simons action. 

Analogously to the description of charged particles interaction by ID integral along Wilson 
loop, we will describe the Casimir interaction of surfaces by corresponding 3D integral. The first 
problem is that for stationary objects the action ([2]) is an integral from t = — oo to t — oo, so 
by analogy to Wilson loop, we should enclose the surface of the integration in t-direction. The 
integration surface for two planes is shown in Fig. 2. This closing procedure can be performed 
both for plane surfaces and for any curved surface in 3-dimensional space. As the result of this 
procedure, so-called Wilson bag [31 H] can be obtained. It can be written as 

e s 

where E is closed 3-dimensional surface in 4-dimensional space-time. The final conclusion is 
that Wilson bag is (by analogy to Wilson loop) observable quantity which gives us Casimir 
energy of the objects, defined by the surface of integration. For two planes we will calculate 
the following object: 

W Bag (R,T) = e us W, (8) 

where 

T 

dxdydz(5(z — R) — 5(z))e 3upa A lJ F pa + 



S(R,T) = J dt 
o 

R 

+ dz 





dxdydt{5{t — T) — 5{t))e Ql/ptT A v . 



o 

And in the limit T — > oo: 



(W Bag (R,T))^Ce- E -^ T . 

The second problem is to rewrite of the Wilson bag in terms of lattice objects (links and 
plaquettes). Using the expansion of link variable: 

Ui, x , p w 1 + igA„{x)a - \g 2 a 2 Al + 0(g 3 ) 
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Figure 3: E cas {R) at /3 = 8.0 and A = 0.0005. 

(here a is the step of the lattice) and plaquette variable: 

CW, « 1 + igF^(x)a 2 - -g 2 a 4 F 2 v + 0(g 3 ), 
we can construct the following expression: 

(1 - Re (U^U^p)) - (1 - RetWO - (1 - Retf I>a . p ). (9) 
It can be easily proved that expansion of (jUJ) by g is: 

^ 3 i^4 p + 0(s 4 ). (10) 

In the main order it gives us one term in the density of Chern-Simons action multiplied by g 2 . 
After this, we replace integration by sum, and multiple this sum by f3 in order to eliminate g 2 
in the last expression. So we obtain a computable expression for Wilson bag. 

Let us consider some results of numerical simulations. You can see on the Fig. 3 the 
dependence of the Casimir energy E cas on R (the distance between planes). 

The parameter j3 is equal to 8.0 and A = 0.0005 for this calculation. The size of the lattice 
is 32 4 . Unfortunately, we can not see here E cas (R) dependence like R~ 3 . All what we can see 
is a distinction of the first point from the other ones which form the plateau. The dynamics in 
the plateau can not be revealed due to numerical errors. Nevertheless, let us make a suggestion 
that this difference of the first point from the plateau is a manifestation of the dependence of 
Casimir energy on the R as R~ 3 (this suggestion follows obviously from dimensional reasons). 
So we have done a fitting of this point sequence by the function 

E C as — G\ + 

R 6 

The values of coefficients C\ and C2 are given on the Fig. 3. Our aim is, of course, the coefficient 
C<i- C\ is an analog of something like self energy of charged particles. Such self energy can be 
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Figure 4: Dependence of C2 on the lattice size. 




Figure 5: Dependence of C 2 on A. 



obtained in Wilson loop calculation. We have performed the same calculations for different 
values of f3. The situation is qualitatively the same. 

Now let us consider the physical meaning of the coefficient C^- A full Casimir energy of 
the interaction between two planes have been obtained in these calculations. This energy is 
expressed in units a -1 , where a is the lattice step. If we take into account that the area of 
the plane is equal to (aN) 2 , where N = 32 is the size of the lattice, we can write the physical 
(dimensional) value of the Casimir energy density: 



E, 



1 Co 1 



CoN~ 



Cas.phys 



a R 3 (aN) 2 (Ra) 



After comparison between this formula and (jSJ) we can conclude that C^A^ -2 should be equal 
to coefficient before R~ 3 in (JSJ) or in (jEJ) (A is small in our calculations by reasons discussed 
later). And this coincidence is rather good. For example, for the calculation with (3 = 8 the 
expression C2N" 2 = —3.27 x 10~ 9 ± 1.1 x 10~ 10 and the theoretical value of this quantity for 
A = 0.0005 is -3.17 x 10" 9 . 

The next step is the analysis of C2 dependences on the size of the lattice A" and on the 
parameter A. The first one is presented on the Fig. 4. We can see here that C2 is proportional 
to A^ 2 with high accuracy. So the calculated Casimir energy of two planes is really proportional 
to theirs area. 

The dependence C%{X) is shown on Fig. 5. In accordance with (jSJ), the calculated coefficient 
at R~ 3 is proportional to A 2 for small A. Unfortunately, present calculation methods can not 
give values of C2 for large A and we have in fact only small A limit. This problem appears 
because the number in the exponent (jHJ) becomes too big and the errors arise for calculations 
with large A. We are studying this problem and hope that some methods of solution can be 
suggested. 
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Figure 6: Dependence of C 2 on (3. 

The last point is the consideration of the continuous limit of our calculations. There are 
two aspects of this problem. First, when we calculate Wilson bag, in fact, only the special 
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Figure 7: Dependence of Ci on (3. 

quantity being approximately equal to Chern-Simons action is calculated (jUJ). Expansions of all 
quantities have been done by powers of g. So the smaller is g, the closer is our lattice Wilson 
bag quantity (jUJ) to its first term in expansion (II Op . Thereby calculations in the limit g> — ► 
or (3 — > oo are the first stage in continuous limit. This stage is shown in Fig. 6 for C<i and 
in Fig. 7 for C\. Additional line at Fig. 6 is the theoretical value of Ci for A = 0.0005. As 
we can see in the limit (3 — > oo the coefficient C2 becomes stable around the right analytical 
value. This stability can be considered as an additional argument for a suggestion that C2 has 
the real physical meaning. As for C\, it hasn't any continuous limit due to permanent growth 
with increasing of f3. 

The second stage is "infinite lattice volume" limit N — > 00. But our method is not very 
sensitive to the lattice size for the plane form of the boundaries. Even for rather small lattices 
(see Fig. 4) the value of the coefficient before R~ 3 is close to the right quantity. There is 
one peculiarity in this continuous limit procedure and we think that this moment should be 
emphasized. The lattice step a disappears from our considerations because the final result 
of our calculation is dimensionless coefficient at -R~ 3 in expression for Casimir energy per unit 
area. So the value of a is not important for us and instead of limit a — > we have limit N — > 00. 

Conclusions 

In this paper we have proposed the numerical method for the Casimir energy calculation based 
on the lattice simulations of QED. We have combined two ideas: the generation of the bound- 
ary conditions by means of the additional Chern-Simons boundary action and the lattice "Wil- 
son bag" concept (the lattice presentation of the closed 3-dimensions surface in Euclidian 4- 
dimensions space). This combination is in fact a lattice definition of the respective quantum 
observable for the Casimir energy and for the Casimir interaction between surfaces. We have 
tested our method in the simplest case of the Casimir interaction between two plane surfaces 
and have achieved a good agreement with the analytical results for this problem. 
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A big advantage of our method is its universality. This method can be applied for following 
calculations: 

• Investigation of the Casimir interaction between surfaces of complicated shapes which are 
interesting from the experimental point of view. 

• Calculation of thermal corrections to the Casimir force. (Unfortunately, this very im- 
portant problem is very difficult for any analytical investigations in physically interesting 
cases.) 

• Study of Casimir effects for non-abelian gauge fields and for fermions which is very essen- 
tial for phenomenological models in heavy-ion collisions and for the models of low-energy 
hadron states [L2l fl~3] . 

These problems are under consideration now. 
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